Effect of particle geometry on phase transitions in two-dimensional liquid crystals 
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Using a version of density-functional theory which combines Onsager approximation and 
fundamental-measure theory for spatially nonuniform phases, we have studied the phase diagram of 
freely rotating hard rectangles and hard discorectangles. We find profound difi^erences in the phase 
behavior of these models, which can be attributed to their different packing properties. Interest- 
ingly, bimodal orientational distribution functions are found in the nematic phase of hard rectangles, 
which cause a certain degree of biaxial order, albeit metastable with respect to spatially ordered 
phases. This feature is absent in discorectangles, which always show unimodal behavior. This result 
may be relevant in the light of recent experimental results which have confirmed the existence of 
biaxial phases. We expect that some perturbation of the particle shapes (either a certain degree 
of polydispersity or even bimodal dispersity in the aspect ratios) may actually destabilize spatially 
ordered phases thereby stabilizing the biaxial phase. 
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I. INTRODUCTION 



Biaxial liquid-crystalline phases were hypothesized in 
1970 by Freiseri Since then many theoretical models 
have described the conditions that promote biaxial or- 
der. However, the experimental search for biaxial phases 
has been unfruitful until very recently, when Madsen 
et ali^ and Acharya et ali^ presented convincing evi- 
dence for the existence of biaxial order in a fluid com- 
posed of boomerang-shaped molecules. In previous work 
Schlacken et al<^ analyzed the effect of particle shape on 
the orientational properties of two-dimensional fluids us- 
ing scaled-particle theory (SPT). They considered hard 
ellipses and hard rectangles (HR) and, for the latter, lo- 
cated a spinodal which corresponds to a tetratic phase, 
where the fluid has two nematic directors oriented per- 
pendicular to each other. Simulation works have been 
performed by Cuesta and Frenkel^ on hard ellipses and 
by Bates and FrenkelS who focused on hard discorectan- 
gles (HDR) and the character of the orientational transi- 
tion from the isotropic phase. Also Lagomarsino et al.- 
found similar behavior in a fluid of hard spherocylinders 
confined between two parallel walls which, in the limit 
of small wall separation, can be considered as a realiza- 
tion of a two-dimensional nematic. These simulations 
have shown that the character of orientational transi- 
tions depends strongly on the particle shape: whereas for 
hard discorectangles the isotropic-nematic transition is of 
the Kosterlitz-Thouless (KT) type^ in the whole aspect 
ratio regime where the nematic phase is stable against 



the solid phase® , for hard ellipses a tricritical point sepa- 
rates regions of continuous from first order orientational 
transition^. 

The density functional (DF) formalism has been 
proven to be a powerful tool in the study of bulk 
nonuniform phases of liquid crystals in three dimen- 
sions, specifically for the model of freely rotating hard 
spherocylinders^"^'*. To our knowledge, in two dimen- 
sions the DF formalism has been applied only to the 
study of uniform liquid-crystalline phases. In the present 
work we will not focus on questions related to the KT na- 
ture of the isotropic-nematic phase transition, but rather 
we investigate by DF theory the occurrence of biaxial 
nematic order in fluids and how this feature depends on 
the geometry of the particles; we consider both HR and 
HDR as interaction models. Also we apply DF theory to 
the study of spinodal instabilities to nonuniform phases 
in order to analyze the relative stability of these biaxial 
phases against spatially ordered phases. 

The results from the theory allow us to conclude that 
a certain degree of biaxiality is present in the case of 
HR. This is truly remarkable, since this order arises in 
a one-component system of particles with uniaxial sym- 
metry. However the bifurcation analysis using the DF 
theory allow us to conclude that the tetratic phase is 
metastable w.r.t. the solid phase. The remainder of the 
paper is organized as follows: In Section II we introduce 
the theoretical models for both the uniform phases (SPT 
in Section II. A) and for the non-uniform ones (DF in Sec- 
tion II. B). Subsection II. C is devoted to the bifurcation 



analysis. Results are shown in Section III while some 
conclusions are drawn in Section IV. 

II. MODEL 

The DF approximation proposed is based on SPT, 
which has proved to be a generic theory which is eas- 
ily extended to any convex particle geometrjii5ii&. Two 
extensions, suitable for nonuniform phases, are proposed. 
Both tend to the SPT in the limit of uniform phases and 
take into account ideas from Fundamental Measure the- 
ory (FMT). This theory has been worked out for a sys- 
tem of hard parallelepipeds in the restricted orientations 
approximation (the so-called Zwanzig model^^) and the 
corresponding phase diagram, including all nonuniform 
phases, was recently calculatedi^. One of the proposed 
extensions recovers the Onsager second-order virial the- 
ory in the low-density limit. The other version is built 
so as to recover the FMT for parallel hard rectangles, 
i.e. restricting the particle orientations to point along 
the nematic director. The construction of a functional 
which interpolates between the Onsager limit and the 
high density one using FMT for hard spheres has been 
proposed by Cinacchi and Schmid for a system of hard 
spherocylinders in three dimensions--. Nevertheless their 
proposal seems to be numerically intractable for nonuni- 
form phases. Here we are interested in constructing a 
workable density functional. 

A. Uniform phases: Scaled particle theory 

The radically different orientational ordering proper- 
ties shown by the two particle models, HDR and HR (see 
Fig. can be understood in the framework of SPT. A 
brief discussion of this theory, as applied to our model 
systems, follows. 

The key quantity to derive the theory is the excluded 
volume between fluid particles of length L (HR) or i + 
(HDR) and width cr, and a test scaled particle of dimen- 
sion Ls and cr^ with fixed orientation Ui((/)i), where <j)i 
is the angle between the uniaxial axis and the nematic 
director. The excluded volume is to be averaged over all 
possible orientations of the fluid particles, 

J d(f>2h{^2)V,l,i{L,,as,<f>12), (1) 

where /i((/>) is the one-particle orientational distribution 
function, 0i2 = 0i — 02 , and V^^^ is the excluded volume 
between two particles which, for rectangles, has the form 

V™i{Ls,as,(f>i2) = {LL, + aas)\sin^i2\+CLCa,{2) 

Cx = VA2 + A2 + 2AA,|cos(/)i2| (3) 

with A ~ {L, a}, whereas for discorectangles 

F™«(L„(T„0i2) = iL + Ls)i<J + as) + ji<J + asf 
+ LL,|sin0i2|. (4) 
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FIG. 1: A system of (a) hard discorectangles in a nematic 
phase, and (b) hard rectangles in a tetratic phase. 

The reversible work required to insert the scaled parti- 
cle with fixed orientation (which coincides with the ex- 
cess chemical potential) is, in the limit of small sizes2£ 
{Ls,(J,) < (L,cr): 

/?Mcxc(0i) ^ /^'^°^(is,a-s,0i) 

= - ln[l-p(14xci)(i.,fT3,0i)], (5) 

where p is the fluid density. In the opposite limit of large 
sizes (Lg, CTg) 3> (i, cr), this work coincides with the ther- 
modynamic reversible work necessary to open a cavity of 
volume equal to the scaled-particle volume (us), which is 
equal to Pv^, P being the fluid pressure. SPT interpo- 
lates between these two limits using a Taylor expansion 
around {Ls,as) = (0,0), where the second-order term is 
Pvg and, finally, the sizes are taken to be those of the 
fluid particles, which results in 

f3fie.M - -ln(l-r;)-f-^ / dcl)2h{^2)V^"H'l>i2) 

+ f3Pv, (6) 

where 77 — pv is the packing fraction and y'-*'-' = 
K;xci(i, cr, 012) — 2u, with V the particle volume. The 
excess chemical potential of the fluid is the angular aver- 
age 

^Mcxc = y oJ0ift.(0i)[/3/icxc(0i)] = -ln(l - ?7) 

+ -^S^+pPv, (7) 
1-77 
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where we have defined Sq = ((V^°^))/(2i;) and ((• • •)) is 
the double angular average. Using the definitions of the 
excess chemical potential and pressure in terms of the 
excess free energy density, ^cxc, a differential equation 
for ipcxc — ^cxcV is obtained, with solution 



-ln(l - 77) + 



1] 



1—7/ 



-Sn 



and, correspondingly, the pressure is given by 



(3Pv = 



(8) 



(9) 



which defines the thermodynamics of the fluid. 

In order to proceed we have to specify a functional 
form for the orientational distribution function, which 
is taken to be a truncated Fourier expansion. In view 
of the symmetry h{4>) = hijr — 0), only even harmonic 
terms have to be included: 



1 



hk cos(2/c</)) 



(10) 



nematic phase, N„, {hi ^ 0) when the following condition 
holds: 1 — (2/37r)?/K^ = 0. Another possibility is that the 
isotropic phase gets unstable when h\ — ^ ^ and 

1 — (2/157r)2/K^ = 0. This bifurcation corresponds to 
a phase with T)^h symmetry, which we call tetratic, Nj. 
The above conditions occur for the following values of 
the packing fractions: 



^N. = (l + Y^4) (15) 



The equality of these two values gives a critical aspect 
ratio K* = (3 + V5)/2 such that, when k < k*, the Nt 
phase preempts the N,i phase. 

In order to calculate the complete phase diagram for 
uniform phases we minimized the full free energy per par- 
ticle for each k with respect to the n coefficients hk- A 
common-tangent construction was used to locate coex- 
istence boundaries in the case of first-order transitions. 
To characterize the orientational order we have used the 
following order parameters: 



fc>i 



h{(j))cos{2k(l))d(p, fc = l,2. 



(16) 



Note that d(f)h{(f)) ~ 1. The cutoff index n is chosen to 
guarantee a small enough value \hn\ < lO"'^. With this 
choice the double angular average Sq for hard rectangles 
becomes 



TT V 



where 



-1/2 



1 

k>l 



(11) 



9k = 



K-1 



K+1 

= (4fc2-l 



\ k = 2j. (12) 



B. Density- functional theory for nonuniform 
phases 

Two different approximations are being proposed for 
the excess part of free energy density <I'cxc(r) which de- 
fines the interaction part of the free energy functional as 
P^cxc — / rfr<I>oxc(r)- The first one will later be used to 
calculate the spinodal instabilities of the HR fiuid against 
nonuniform phases. 

As an important ingredient of any density functional 
approximation for freely rotating anisotropic particles, 
we need to impose the Onsager form for the excess part 
of the free-energy density in the low density limit: 



We have defined k = L/a (the aspect ratio). The equi- 
librium orientational structure of the fluid is obtained by 
minimizing the free energy per particle, whose ideal part 
reads 



t£id 



hu] — 1 



d(j)h{(l)) In [7r/i(0)] 



(13) 



whereas the excess part (pexc/?7 is obtained from Eqn. 
A bifurcation analysis for small values of the hk coeffi- 
cients leads to the following expression for the difference 
between free energies A(/3 = (/Sjv — of nematic (N) and 
isotropic (I) phases: 



'7 



hi 



IStt 



hi 



(14) 



where k± = k^/^±k^^^^ and y — rj/ {l — rj). Note that the 
first term is dominant with respect to the second. The 
isotropic phase gets unstable with respect to the uniaxial 
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X /(ri,r2,( 



.)■ 



c^</>2p(ri,^i)p(r2,02) 
(17) 



where /(ri, r2, 0i, ^2) is the Mayer function between two 
particles with orientations 0i and 02- Another impor- 
tant feature is that the functional capture correctly the 
high-density limit. For this purpose we will use the FMT 
to include the remaining part of the density dependence 
and ensure, in this way, the high density limit. This de- 
pendence enters through the weighted densities ^^(r)^^. 
Finally we impose that the scaled particle theory be re- 
covered in the uniform limit. 

As a first proposal we write the following form for the 
excess part of the free energy density: 

$...(r) = -no(r)ln[l-n2(r)]-^^^^(^ 

1 — n2(r) 



(1 



(r)<i>Of(r), (18) 
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where the n„'s are defined 



(a) 

P* OJy 



(r), i.e. 

they are convolutions of the averaged density profile of 
rectangles, p(r) = J d(j)p{r, (p), with some weighting func- 
tions given by the expressions 



,(0) 



,(2) 



(r) = e 



e 



\v\ 



(19) 
(20) 



where 5{x) and Q(x) are the Dirac delta function and the 
Heaviside function, respectively. An important property 
of the above proposal is that it recovers the SPT theory 
in the uniform limit. Note that weighting functions are 
the characteristic functions which define the support of 
the particle geometry (a rectangle parallel to the y axis) 
and its total surface areaii (their integrals are two of 
the fundamental measures of the particles). Note also 
that, in order to numerically implement the theory, the 
simplest choice is to take a system of parallel HR as a 
reference system. Another possibility is to use as a ref- 
erence system a fiuid of parallel hard squares with par- 
ticles having the same volume as that of the rectangles. 
This choice may be more adequate to describe isotropic 
nonuniform phases such as the plastic solid (for which 
the lattice period is the same in the x and y directions), 
whereas the former choice should be better for the ori- 
ented solid phases. Both reference systems will later be 
used to calculate the spinodal instabilities to nonuniform 
phases. 

A possible improvement of the theory would be to get 
rid of the reference system of parallel particles and gener- 
alize the definitions of weighting functions H19|) and H20|) 
to describe freely rotating particles: 



4 V2 



L 



c.(2)(r,0) = e(|-|x^|)e[^-|r^ 



(21) 
(22) 



where = a; cos ^-fy sin and = ycos(j)—xshi(j) are 
the new Cartesian coordinates in a system rotated by the 
local angle </> of the particle, measured w.r.t. the nematic 
director (the y axis). Thus the weighted densities n^'s 
are now calculated from 

n„(r) = j dv' j dcjy'piv + r', 0')<nr', <P')- (23) 

The new expression for <i>oxc is now given by Eqn. IjlSj) 
where the last term is substituted by 



d4)i 



(l-n2)-i*<) (ri,0i)p(ri,0i) 



j dY2 j d(/)2p(r2,02)/(ri,r2,0i,(/)2). 



(24) 



This term has to be modified due to the fact that the term 
(1 — 712)"^ is now convoluted with the weight u}^\v, (f>), 



which explicitly depends on (f>. This approach will not be 
pursued any further and we leave it for future work. 

As can be seen the prefactor of the function $2cc^ 
Eq. H18|) depends on the density profile defined in a rect- 
angular region centered at r with dimensions (2i) x {2a). 
This is not the usual case in the formulation of the FMT 
from first principles, but it is necessary because the term 
^^^^ depends on the Mayer functions which have the 
range of the excluded surface area. 

In order to test the performance of this version of DF 
we will apply it to the study of freezing in a parallel 
hard-square (PHS) fiuid and compare the results with 
the FMT derived inil. This test is motivated by the fact 
that the close packing limit is easier to reach in this sys- 
tem compared with the freely rotating case and thus the 
possible divergences inherent to the particular functional 
structure can be more easily elucidated. 

The density profile of the solid phase was parameter- 
ized using Gaussian peaks centered at the sites of a square 
lattice and the fraction of vacancies was included through 
a normalization factor in the profile. The minimization of 
the functional was carried out with respect to the Gaus- 
sian width and occupancy probability v (one minus the 
fraction of vacancies). The results are shown in Fig. [21 
where the fiuid and solid equations of state (Fig. 12^) and 
the occupancy probability vs packing fraction (Fig. I^b) 
are plotted. Comparison is made between results from 
the present theory and FMT. The fluid branch is the 
same because, as already pointed out, our theory recov- 
ers SPT. Note that the new version stabilizes the solid 
phase at a lower packing fraction rj w 0.50 with a frac- 
tion of vacancies 1 — v = 0.085, to be compared with 
the FMT result (r/ w 0.53 and 1 - 1/ = 0.154). The lat- 
tice parameter, d, of the solid phase at the bifurcation 
point, as calculated from the mean density p = I'/d'^, is 
equal to 1.353 from the present version and 1.254 from 
FMT. All of these results indicate that both function- 
als produce similar results and, what is more important, 
the new functional is regular for this kind of density pro- 
files (note that a free minimization would be required to 
completely settle this point). 

Another possibility to construct a DF for freely rotat- 
ing hard rectangles is to impose that this functional re- 
cover the Fundamental Measure Functional (FMF) when 
the density profile is chosen to have the form p{r, (p) — 
Py{r)5{(p), i.e. the system of parallel hard rectangles. 
Note that, in this limit, we should recover the equation 
of state of the solid corresponding to a system of parallel 
hard squares, shown in Fig. |21i. A free-energy density 
that conforms to this criterion is proposed in the Ap- 
pendix. 



C. Bifurcation analysis of nonuniform phases 

At fixed chemical potential fiQ corresponding to a ne- 
matic phase characterized by the one particle distribution 
function pQh{(j)), the minimization of the grand potential 
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FIG. 2: (a) Equations of state for PHS's as obtained from 
FMT and the proposed version of DP. Dotted line: fluid 
branch; dashed line: solid branch from FMT; and solid line: 
solid branch from the proposed version. The full square and 
the circle indicate the corresponding bifurcation points, (b) 
The occupancy probability vs packing fraction along the solid 
branches. Dashed line: FMT, and solid line: present theory. 



leads to 



where the hats over the functions ^ and C indicate 
Fourier transforms. 

Inserting the truncated Fourier expansion of the func- 
tion |(q, 0), 



|(q,0) =5;]6(q)cos(2fc0) 



(28) 



A;=0 



in Eqn. (|27|l . multiplying the latter by cos{2j(f) and inte- 
grating in (j) from to 27r we obtain the following algebraic 
equation: 



(/-Por)u = 0, 



(29) 



where / is the identity matrix, T is a (n -|- 1) x (n + 1) 
matrix with elements 



1 



d4>h{(t)) 



7r(l + 5jo) 

X / d<?!)'cos(2j(/))cos(2fc0')C'(q, (/),(/)'), (30) 



and u — [Co(q), ■ • • , ?n(q)]"^ is a n -I- 1 vector. 

The lower value of po and the vector q for which there 
exists a non trivial solution to Eqn. (|29() can be calcu- 
lated from 



i7(po,q) = det(/-poT^) = 0, (31) 
Vi/(po,q) = (32) 

Both equations will be used later to calculate the spin- 
odal instabilities from the isotropic or nematic phases to 
possible nonuniform phases. 



0) = Poh{<t>) exp 



-/3 



6p{r, (j)) 



Mo 



(25) 



Assuming a density profile of the form p{r, (f) = 
p^hijp) -f$(r, (f) for the incipient nonuniform phase which 
bifurcates from the nematic phase, a Taylor expansion of 
Eq. (|25l) up to first order in ^ gives 



P(i", <t>) = Po^('/>) 



1 



j dv' j d0'C(r-r',0,0')C(r',0') 



, (26) 



where C(r — r',(j>,<j)') is the direct correlation function 
of the nematic fluid calculated as minus the second func- 
tional derivative of PTc-mc with respect to the density pro- 
file evaluated at ^ = 0. Eqn. (|26|) can be rewritten in 
Fourier space as 



|(q, 0) - pM) / d</)'C(q, 0, (A')l(q, 4^') = 0. (27) 



III. RESULTS 

A. Uniform phases 

In Fig. 13 we show the phase diagram for hard rectan- 
gles with 1 < K < 25 as obtained by applying the SPT 
formalism. For 2.21 < k < 5.44 the I— N^i transition 
is of first order. For n < 2.21 a tetratic nematic phase 
Nt begins to be stable; its stability region is bounded 
below by a second-order I— Nt transition. The region 
is bounded above by a Nt— N,i first-order transition for 
1.94 < K < 2.21 and by a second-order transition for 
1 < K < 1.94. This means that there exist two tricritical 
points: one at k = 5.44, already predicted by Schlacken 
et al."', and a second one at k = 1.94, which can only be 
predicted by a proper bifurcation analysis from the Nj to 
the Nu phases. Fig. ^ shows the two order parameters 
[Eqns. IjlGIl ] along the transition lines to orientationally 
ordered phases, for the systems where these parameters 
are different from zero. Orientational distribution func- 
tions for the Nu and Nt phases at coexistence for k = 2 
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are shown in Fig. [S| Note that the Nt phase presents two 
peaks of the same height (consistent with this phase be- 
ing invariant with respect to the D^h symmetry group). 
By contrast, the N„ phase also presents a secondary peak 
(indicating a certain degree of biaxiahty in this phase) , 
but of a smaUcr height. If k is increased the height of 
this secondary peak reduces compared to the main peak, 
eventually disappearing. 
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FIG. 3: Phase diagram for hard rectangles in the rj — k plane. 
The K axis is represented using a logarithmic scale. Dashed 
lines indicate second-order transitions, whereas full lines are 
binodals. The shaded region is the two-phase region where 
phase separation occurs. Labels indicate the three phases 
involved: isotropic (I), uniaxial nematic (N„) and tetratic 
nematic (Nt). 




FIG. 4: Order parameters qi and q2 along the transition lines 
to orientationally ordered phases. The full line corresponds 
to the qi parameter, whereas the dashed lines correspond to 
the q2 parameter. 

We now show the results for hard discorectangles. We 
have implemented the same calculational scheme and ob- 
tained the phase diagram. Only two phases are involved: 
the isotropic phase I and the uniaxial nematic N phase. 



The phase boundary r] — k between the two turns out to 
be of second order always^, and is given by the equation 



2 (k-1)/tt 



(33) 



4(«-l), 



where k = {L + a)/a and now L is the length of the 
rectangular part of the particle and a is its width. 

A key difference between hard rectangles and hard dis- 
corectangles is that in the former case the excluded vol- 
ume expression contains a cosine term which is absent in 
the latter case. This can be seen from Eqns. and 
(|H) : substituting Lg = L and Cs = f we obtain 



K^«i(L,a,0i2) 



(i^ +cr-J|sm0i2| 

2Lcr(l-h|cos(?!)i2|) (34) 



whereas for discorectangles 



HDR 



(i, Cr, 012) 4L(T -I- TTcr^ + ^2 



(35) 
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FIG. 5: Angular distribution functions of hard rectangles of 
the Nu (solid line) and Nt (dashed line) phases at coexistence 
and for k = 2. The dotted line represents the distribution 
function of hard discorectangles for the same k and for r) = 
0.91. 

This basic difference translates into the gk prefactors 
in Eqn. Ijlll) being different according to the parity of 
the index for hard rectangles [see Eqn. dJ], whereas 
for hard discorectangles these coefficients have the same 
form. This feature again shows up in the orientational 
distribution function which, as pointed out above, ex- 
hibits a secondary maximum at = 7r/2 in the case of 
rectangles, a peak which is absent in the case of discorect- 
angles, see Fig. |S1 In this figure we also show the angu- 
lar distribution function corresponding to the N phase of 
hard discorectangles for the same k (and for 77 = 0.91). 

In order to compare the equations of state of hard 
rectangles and discorectangles for the same aspect ra- 
tio K = 5, the reduced pressure vs packing fraction is 
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shown in Fig. |B1 It can be seen that the two equations 
of state show a similar behavior, except in the neighbor- 
hood of the corresponding phase transitions, which are of 
first and second order for rectangles and discorectangles, 
respectively. 




FIG. 6: Equations of state of hard rectangles ( solid line) and 
discorectangles (dotted line) for k — 5. The arrow indicates 
the location of the second order phase transition. 



B. Nonuniform phases 

To elucidate the absolute stability of the tetratic ne- 
matic phase (see Fig. we have carried out a bifurca- 
tion analysis following the lines described in section lTl CI 
i.e. solving Eqns. and (|32|l in order to find the val- 
ues of the packing fraction and wave vector for which the 
uniform phases destabilize to nonuniform phases. The 
results are shown in Fig. [7| where we plot all spinodal 
lines of the TNu, the I-S (S being a spatially ordered 
phase, whose structure will require further analysis) and 
the Nu-S transitions. Note that Fig. [Tjis in fact the same 
as Fig. 121 but spinodal lines from the above bifurcation 
analysis have been superimposed and, at the same time, 
transition lines in Fig. |31which turn out to be preempted 
by nonuniform phases have been dropped. 

As can be seen in Fig. [3 for k < 5.5 the I phase desta- 
bilizes directly to the S phase, so the tetratic nematic 
phase is not stable. Nevertheless, a tetratic order is pos- 
sible in spatially ordered phases. Taking into account 
that the l-S transition is probably of first order a com- 
plete functional minimization with respect to two dimen- 
sional density profiles is required to settle this question. 
In Fig. |7|the nonuniform phase spinodals are calculated 
using two reference systems: parallel rectangles and par- 
allel squares. As was discussed in section III Bl the last 
one should better describe the I-S transition, whereas the 
former should be more appropriate for the Nu-S transi- 
tion. In other words, the dotted line should be more 
reliable for k < 5.5, whereas the solid line is expected 
to be more accurate for k > 5.5 (note that both refer- 



ence systems give similar results in the first region and 
naturally these results coincide exactly at k = 1). 

If we use the density functional for freely rotating par- 
ticles without a reference system of parallel particles (see 
section lTl B|) . the fluid-solid transition is expected to shift 
to higher densities. The reason for this behavior is that a 
reference system composed of parallel particles enhances 
a crystalline local order which is absent in the isotropic 
phase. Nevertheless, we expect that, within the present 
version of the theory, the tetratic phase still remains un- 
stable w.r.t. the solid phase. Calculations using this 
functional is a task in progress and will allow us to fi- 
nally settle this question. 
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FIG. 7: Phase diagram of hard rectangles including the spin- 
odal lines of the I-S and Nu-S transitions. The solid and 
dotted lines represent the spinodal curves calculated from the 
functional obtained from Eqn. 11811 . using a reference system 
of parallel rectangles and parallel squares, respectively. For 
the latter case the volume of the squares were taken to be the 
same as that of rectangles. 



IV. CONCLUSIONS 

The main conclusion of the present work is that dif- 
ferent particle geometries may account for the differ- 
ent nature of the orientational phase transitions involved 
(continuous versus first order with a tricriticcal point). 
But, even more important, particles possessing uniax- 
ial symmetry (hard rectangles) may lead to a phase 
with a dramatically different symmetry: the tetratic 
phase. This finding is based upon minimization of a 
density-functional theory constructed using ideas from 
SPT. However, the tetratic phase has been found to be 
only metastable with respect to a spatially ordered phase. 
This result was obtained using a bifurcation analysis of 
a version of density-functional theory proposed in the 
present paper. The theory combines the Onsager func- 
tional with the FMF for hard rectangles. A word of 
caution is in order here: the proposed functional does 
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not reproduce the exact third virial coefficient, which 
is expected to be of vital importance at high density, 
and these neglected correlations should play a role of 
paramount importance in the stabilization of a phase 
with tetratic symmetry. It might well be that a full the- 
ory, incorporating correlations between three particles, 
could stabilize the tetratic phase even further and predict 
it to be stable at densities below that at which crystal- 
lization occurs. Nevertheless, it is interesting to note that 
the present theory already confirms that an oriented fluid 
of hard rectangles presents some degree of biaxial order 
for small aspect ratios. The theory could also be used 
to predict the stability of a number of two-dimensional 
phases that occur in Langmuir monolayers'^. 

Previous work on polydisperse mixtures of hard 
spherocylinders^2i24, or parallelepipedsSS. in three dimen- 
sions, have shown that polydispersity enhances the sta- 
bility of the uniform phases w.r.t. nonuniform phases. 
While columnar order admits a higher degree of poly- 
dispersity, the smectic or solid phases normally have a 
terminal polydispersity beyond which these phases are 
no longer stable^?,. This fact can be easily understood 
if we take into account that: i) any periodic packing is 
difhcult to attain in a mixture of particles with different 
characteristic lengths and ii) mixing increases the ideal 
mixing entropy of uniform phases and, as a consequence, 
decreases the Helmholtz free energy. Thus, some degree 
of polydispersity (either unimodal or bimodal) will cer- 
tainly increase the density of the transition from the fluid 
to a spatially ordered phase in a polydisperse mixture of 
hard rectangles. Is this effect enough to stabilize the 
tetratic phase? This is an open question that we will try 
to elucidate in a future work. 
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the uniform limit. Our proposal for the excess free-energy 
density is 



$cxc(r) = -no{r) ln[l - n2(r)] 



'2E 

fj,=x,y 



(r) 



1 - n2(r) 



/ fir'/o(r')/ref (r - r' 



Ty (36) 



where /icf(i") — —^{(^ ~ \x\)Q{L — \y\) is the Mayer func- 
tion of two parallel hard rectangles, <I>Jj^c is the excess 
part of the Helmholtz free-energy density obtained from 
the FMT for parallel hard rectanglesii, i.e. 



rcf 



-raoln(l - 712) + 



1 - 712 ' 



(37) 



while the new weights are 
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(38) 
(39) 



Note that the uniform limit of $0x0 (r)w [see Eqn. H36|l ] 
coincides with the expression (jSJ obtained from the SPT. 

Now we will show that the free-energy functional ob- 
tained from Eqn. (|36|) reduces to the FMT for the case 
of parallel rectangles. In this limit it is easily shown that 



$c".^^(r) 



-Ar) 



/ dr'p(r')/ref (r - r') 2' 
The above limit, together with the fact that 



(40) 



dr 



1 — 712 



(r)p(r) 



1 - 772 (r) 



(where n ^ v) allow us to conclude that 



Appendix 

In this Appendix we present an alternative version of 
DF which recovers the FMT for parallel HR and SPT in 



dr$oxc(r) 



dr$-f(r). 



(42) 
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